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RELATIVISTIC MODELS OF THE UNIVERSE 
WITH PRESSURE EQUAL TO ZERO AND 
TIME-DEPENDENT UNIFORMITY 

By Windsor L. Sherman and Sylvia A. Wallace 
Langley Research Center 

SUMMARY 

The previous work on the zero-density and approximate models of the universe, 
published in NASA TN D-2601 and TN D-3047, is summarized. The zero-density model 
is improved so that it is a useful approximation of the exact model to greater values of 
the red shift, and regions of applicability of the approximate model are determined. 
Regions where closed-form solutions of the exact model exist, other than zero-density 
and zero- cosmical- constant models, are defined and the solutions are given. The results 
of calculations are presented that indicate, for the data considered, all the models, sim- 
plified and exact, appear to fit the data equally well even though there is wide variation 
in the density and acceleration parameters. In addition, the results of the calculations 
indicate that density and acceleration parameters should be determined as a pair from 
observational data. It is concluded that data consistent with that for galaxies at red 
shifts greater than one-half are required before a model universe can be selected and a 
simplified model defined. 


INTRODUCTION 

One of the uses of the theory of model universes is the analysis of observational 
data to obtain an indication of the structure of the universe and how the universe has 
evolved. There are several theories of gravitation that provide theories of model uni- 
verses. The discussion in this paper is limited to models of the universe based on 
general relativity. The key assumptions made are that the m'odels have time-dependent 
uniformity and that the pressure is low enough to be assumed negligible when compared 
with the density term. The exact equations for the relativistic models with zero pres- 
sure and time uniformity are derived. The exact equations contain elliptic integrals and 
therefore cannot, in general, be integrated to obtain analytic solutions. The use of the 
simplified models, presented in references 1 and 2, is discussed and additional material 
relative to these models and their use is included. (Appendix A gives the regions of 
usefulness of the approximate model.) In addition to the zero-density and approximate 
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models and those models in which the cosmical constant is zero, there are other models 
that have solutions in terms of simple functions. McVittie (ref. 3) has discussed these 
models from the point of view of special solutions of the elliptic integrals. In this paper, 
an algebraic approach is used that permits the establishment of the combinations of the 
density and acceleration parameters that give additional simple function solutions. 

The differences between the meaningful simplified models and the exact models 
are discussed, and the results of calculations with the various models to determine the 
density and acceleration parameters are compared. Appendix B lists the symbols used 
in the analysis. 


RELATIVISTIC MODEL UNIVERSES 


The space-time metric for relativistic models of the universe that have time- 
dependent uniformity is 


ds 2 = dt 2 




dr 2 + r 2 dg 2 + r 2 sin 2 g d<£ 2 



( 1 ) 


where ds has the dimensions of time along the line element; c is the speed of light 
in vacuo; t is time; r, 0, and 0 are the dimensionless coordinates of a point in the 
metric subspace; k is the curvature constant; and R is the scale factor that describes 
the geometrical history of the universe and has the dimensions of length. The relation 
between the space time of experience, the metric space, and the physical content of space 
time is provided by the field equations of general relativity: 

® l±v ~ 2 ^=^ - 2Ag ]UZ ,j = -kc^T^ u (2) 

where is the contracted Riemann- Christoff el tensor; E, the curvature scalar; 

g„ j,, the metrical tensor; A, the cosmical constant; k is a constant = and G is 
h- c 2 

the universal constant of gravitation. The material energy tensor T^ v is defined as 
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( 3 ) 


for a single stream isotropic fluid. In equation (3), p and p are the density and 
pressure, respectively, and V^ = — The subscripts p and v appearing in equa- 
tions (2) and (3) are the usual tensor indices and can assume values 0 to 3. 


The field equations (eq. (2)) are written with the cosmical constant. Einstein’s 
original formulation was 
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r 


E /xy - ^ E Spv - -KC^T !± v 


( 4 ) 


The eosmical constant was introduced by Einstein in reference 4 in order to obtain a 
stationary model universe which at that time was indicated by observation. After obser- 
vation had established the expansion of the universe and it had been shown that expanding 
universe solutions existed for A = 0 as well as A 4 - 0 (ref. 5), Einstein discontinued 
the use of the eosmical constant on the grounds of logical economy (ref. 6) and maintained 
this view through his last discussion of the question (ref. 7). The question of the cosmi- 
cal constant is not as simple as the foregoing discussion would indicate, and Einstein him- 
self was not consistent in his use or rejection of the eosmical constant. When working 
with the electromagnetic field (ref. 8), Einstein found and used the eosmical constant as 
"a constant of integration," whereas when working with the ponderable matter field, he 
regarded it as a constant peculiar to the theory. Weyl (ref. 9) presents a rigorous math- 
ematical derivation of the field equations and shows that the eosmical constant is a neces- 
sary term in the field equations. McVittie (ref. 10) and Heckmann (ref. 11) both support 
the retention of the eosmical constant on mathematical grounds. Because the mathemat- 
ical evidence indicates that the eosmical constant is a necessary term in the field equa- 
tions and the rejection of this term seriously limits the generality of the theory of rela- 
tivity, the eosmical constant is retained in this study. 


The substitution of equations (1) and (3) into equation (2) and carrying out the 
indicated operations gives the well-known form of the Einstein field equations for an 
isotropic homogeneous universe 


^ + =9 + « P = - A 

R R 2 R 2 


(5a) 


R.2 _ kc^p _ _ kc^ A 

r2 3 r2 3 


(5b) 


where R is the scale factor of the universe and the dots over the symbols indicate dif- 
ferentiation with respect to time. The solutions R(t) of these differential equations 
describe the geometrical history of the universe. There are many possible solutions of 
equations (5), but it is required that the solutions have physical significance. Since 
observation indicates the presence of matter and radiation in the universe, only those 
sets of k, A, and R that give p 5 0 and p ^ 0 are of interest. 

An examination of equations (5) shows that A has the dimensions of sec -2 . The 
eosmical constant can be greater than, less than, or equal to zero. When A is less 
than zero, it represents an acceleration term that acts in conjunction with gravity forces 
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to retard expansion. When A is greater than 0, it represents an acceleration that 
opposes that of gravity in the retardation of the expansion of the universe. 

Like A, k can be greater than, less than, or equal to zero. The sign of k 
determines the type of space that is being considered. For k > 0, the space is spherical; 
for k < 0, the space is hyperbolic; and for k = 0, the space is Euclidean. In the deriva- 
tion of equations (5), k was adjusted so that it can be considered a space- curvature 
constant that can assume a value of 1, -1, or 0. 

In equations (5) p and p are functions of time; thus, equations (5) are not deter- 
minate since three unknowns must be determined from two equations. Solutions to equa- 
tions (5) can be obtained only if some assumption is made with regard to p and p. 
Before such an assumption is made, however, it is of interest to look at the magnitudes 
of the density and pressure terms. The pressure is composed of two terms, one due to 
matter and the other due to radiation. In a low-density fluid the pressure due to mate- 
rial density is given by two-thirds of the energy density due to random motions which 
reduces to p 0 v 0 3 where v 0 is the random radial velocity. The expression for total 

2 aT^ 

pressure can now be written as p = p 0 v 0 + where p 0 is the material density, 
a = 7.564 X 10" 1^, and T is the radiation temperature in degrees Kelvin. The term 
aT^/3 is the radiation pressure. A good value for p 0 , based on observation, is 
10~30 g/cm^. (See refs. 12 and 13.) The random radial velocity has not been accu- 
rately determined but is probably less than 3 X 10^ cm/sec (ref. 14). These values of 
p 0 and v 0 give a pressure due to material density on the order of 10" 15. The radia- 
tion temperature of intergalactic space, for the present epoch, is about 3° K (refs. 15 
and 16) and corresponds to a radiation pressure on the order of 10" 13. These results 
indicate that for the present epoch, radiation pressure dominates the pressure term. 


Density is also composed of two terms, the material density (that is, the density 
due to matter) and the density equivalent of radiation (that is, aTVc^). For a radiation 
temperature of 3° K the density equivalent of radiation is on the order of 10" 33 and thus 
smaller by a factor of 10 3 than the density due to matter. 


From equations (5) the relation between pressure and density is 



+ 3® r3 = ° 


( 6 ) 


This equation can be written as 



The results of the analysis 


of density and pressure for the present epoch, when applied to the coefficient of dR 3 /dt, 
indicate that p/c^ will be less than p by a factor of 10^ ; thus, pressure can be 
assumed to make a negligible contribution to the value of dR3/ dt when compared with 
density. Thus, in the models of the universe discussed herein, the pressure is assumed 
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to be zero. As a consequence of the assumption of zero pressure, integration of equa- 
tion (6) gives pR 3 = Constant = p 0 R 0 ^ and gives density as function of R and 


R 


Of 



The subscript o indicates present values of p and R. 


Under the assumption of zero pressure, equations (5) become 


2R R^ 
R R 2 



R 2 _ kc^p _ _ kc 2 A 
R 2 3 R 2 3 


(7a) 

(7b) 


Equations (7) are the equations that are used in the development of model uni- 
verses. When equations (7) are written for the present epoch and the following defini- 
tions are substituted 



R o 

R 0 


— -q 0 H 0 


2 


a o = 


4irGp o 

3H q 2 


it is found that 


A = 3H O 2 (0 O - q Q ) 
kc 2 2 

zr 2 - H o v 3o o - % - 1 
iV 


( 8 ) 

( 9 ) 


where the density parameter o Q and acceleration parameter q Q are pure numbers and 
the Hubble parameter H 0 has dimensions of sec - *. If ct q , q Q , and H 0 can be deter- 
mined from observation, equations (8) and (9) can be used to determine the type of space 
curvature and the value of the cosmical constant. If the curvature constant and cosmical 
constant are known, the type of model universe can be determined. For zero-pressure 
models, the models that correspond to k = ±1 or 0 for A % 0 can be determined from 
equation (7b) without reference to o Q and q 0 . 

To determine R(t) for the various combinations of k and A, equation (7b) is 
solved for time; the solution yields 
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R 1 / 2 dR 


( 10 ) 


*0 - 




- kc 2 R + 2H Q 2 a 0 R 0 3 


and the inversion of the time solution gives R(t). If the integration indicated in equa- 
tion (10) is carried out, the solution is in terms of elliptic functions. However, study of 
the cubic in equation (10) provides information on the behavior of R in different models 
This type of analysis has been made by both Robertson (ref. 17) and Bondi (ref. 18), and 
table I summarizes the results of these studies for model universes in which p > 0 and 
P = 0. 

The material presented in table I is self-explanatory for k = -1 and 0 and k = 1, 
A = 0 since for a given combination of k and A, only one type of model universe is 
possible. When k = 1 and A > 0, six different models of the universe can occur. The 
specific model depends on whether A is greater than, equal to, or less than a specific 
value called the critical value (designated Ae). If the cubic in equation (10) is set equal 
to zero and solved for A, it is found that 

A = 8kc oo o (n) 

R 2 R 3 

and the critical value is that value of A for which dA/dR = 0 is given by 


A e = 


k 3 c6 


9H 0 4r 0 6 CTo 2 


( 12 ) 


The value of R associated with A e called R e is given by 


R e - 


3H 0 2 R 0 3 °q 

kc2 


(13) 


If A > A e > 0, the model is Le Maitre’s model (indicated as model 9 in table I). When 
A = A e , three different models are possible. The model in this case depends on the 
relationship of R 0 to R e . When R 0 < R e (model 12), a model that expands asymptot- 
ically from R = 0 to R = R e occurs; when R 0 = R e , the Einstein static model 
(model 10) is indicated; and when R 0 > R e , the Eddington-LeMaitre model (model 11) is 
the model that would represent the universe. When A < A e , the cubic has two roots, 

Rj and R 3 , Rj < R 2 ; when R 0 < Rj, then model 13 occurs; and when R 3 = R 0 , the 
model is model 14. 


Some interesting points about the models of the universe are shown in table I. In 
most cases a singular point at R = t = 0 indicates that expansion started from a zero 
radius. This singular point is purely mathematical and the theory of model universes 
does not provide information concerning the physical state of the universe. There is 
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TABLE I.- SUMMARY OF MODEL UNIVERSES 



A <0 


Oscillating 


R 
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TABLE I.- SUMMARY OF MODEL UNIVERSES - Concluded 


Type 


k = 1 


R(t) 




some indication (ref. 3) that the singularity may have been introduced by the assumption 
of zero pressure, and if a finite pressure model is considered, R > 0 at t = 0. The 
other point of interest concerns the models of the universe in which A = 0. When 
A = 0, equation (10) can be integrated by simple functions and the properties of the A = 0 
models are given by simple analytical expressions. However, by equation (8) a Q must 
be equal to q Q where A = 0. This condition is most restrictive; in addition to 
requiring the equality of a 0 and q Q , the universe cannot be accelerating. Accelera- 
tion requires negative q Q and negative q Q gives negative density; the A = 0 models 
and negative density cannot exist in the physical universe. The mathematical simplicity 
of A = 0 is most attractive; however, the restrictive nature of the assumption lends 
support to the mathematical reasons previously cited for retaining A in the field equa- 
tions. If A is actually zero, the analysis of observational data will show this condition 
and firm justification will exist for setting A = 0. 


Connection Between Theory and Observation 

The classes of observational data that can be used to determine a model of the 
universe are apparent magnitude, red shift, angular diameter, and number counts. The 
isophotal diameter is not considered because of difficulties associated with its use. (See 
ref. 14.) In order to determine a model of the universe from observational data, it is 
necessary to write down an equation that expresses a relation of two or more of the 
classes of observational data as a function of a Q and q 0 . From this equation, expres- 
sions relating apparent magnitude and red shift, angular diameter and red shift, and num- 
ber counts and red shift have resulted. These relationships, in their most general form, 
are: 


Between red shift and apparent magnitude (refs. 1 and 3): 

m - K = 5 log R 0 (l + z)S(co) + M 0 + AM 0 - 5 
Between angular diameter and red shift (ref. 19): 

0 = S(L+ z -). 

R 0 S(cu) 

Between number of galaxies and red shift (refs. 19 and 20): 

N(m) = - — Q- -^sinh \[-ku> cosh\/--ka> - \pk(vj 


where 


(14) 

(15) 


(16) 


m apparent magnitude 

K red- shift correction 
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z red shift 

M 0 absolute magnitude of equivalent local source 

AM 0 evolutionary correction to M 0 

<f) angular diameter j 

! [ 

£ linear diameter ! 

N(m) number of like sources brighter than a given magnitude 

n number of like sources per unit volume 

Q number of square degrees in the celestial sphere 


Equation (14) is the red-shift — magnitude relation and the argument of the log term 
R 0 (l + z)S(co) is the luminosity distance D^, the distance that light travels from the 
source to the observer. Equation (15) is the angular- diameter — red- shift relation and 
the denominator of the right-hand side R 0 S(w) is the distance by measuring rod D c . 
Equation (16) is the number count relation and, as indicated in reference 21, it is not too 
useful with presently available data. 

An inspection of equations (14), (15), and (16) shows that R 0 , w, and S(a>) are 
functions of the radial metric variable r of equation (1), and they need to be expressed 
as functions of a Q , q Q , and z in order to make the equations useful for the analysis 
of observational data. Equation (9) gives R 0 as a function of a Q and q Q . The vari- 
able a) and the function S(w) can be obtained from equations (1) and (7) as functions 
of <7 0 , q Q , and z. 

Equation (1) is first rewritten with the aid of the transformation 


R 0 

£ = sinh 


4-kco = 


R 0 r 


1 


+ 


kr 2 

4 


so that the metric (eq. (1)) becomes 
ds 2 = dt 2 - 


CM 

d£ 2 

c 2 R 0 2 

i-*4 

R ^ 


£ 2 (d0 2 + sin 2 0 d$ 2 ) 


(17) 


(18) 


If the coordinates are chosen so that the origin is at the observer, by spherical symmetry 
6 and <fi are constants for any given light ray. The coordinates are so adjusted that 
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6 = (p = 0. Light travels along a null geodesic and ds 2 = 0; thus, equation (18) can now 
be written as 




c^ = 
R 


Rr 


Differentiation of equation (17) gives 


dw = 


K 0 


AL 


Hoi 1 


. Ml 

Rn 2 


If this definition of dw is used, equation (19) can now be written as 


, dt 
da) = c — = 
R 




Roll 1 


. Ml 

c 

W 


Rr 2 


and integrating yields 


From equation (20), 


and 


CO 


= -Lsinlf 1 ^ J-=c f t0 
y-k R o J t 


dt 

R 


S(co) = — sinh \f-kco 

fk 


CO 


= c r to — = r R ° dR 

J t R Jr RR 


(19) 


( 20 ) 


(21a) 


(21b) 


Through the use of equation (7b), the definition of cr n , and the relation between p 

( R 3 \ 

and p Q for zero pressure models (p 0 = p — gj, equation (21b) can be written as 

Rq 


U) 


-c 


dR 


R 1 / 2 ^2H o 2 R o 3 ct 0 - kcR 2 +|R S 


( 22 ) 


Equation (22) can be expressed in the terms of o Q , q Q , and the red shift z through 
the use of equations (8) and (9) and the relation R 0 /R = 1 + z. Equation (22) can now be 
written 

dz (23) 


CO 


=^r 

H 0 R 0 


° [ 2cr o z3 + ( 3cr o + % + !) z2 + 2 (% + 1 ) z + ^ 


1/2 
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Equations (14) to (16) can now be expressed in terms of ct q , q Q , H 0 , and z by the use 
of equations (9), (21a), and (23). 

Integrations To Obtain on 

Equation (23) is an elliptic integral and cannot, in general, be integrated in terms 
of simple functions. There are certain assumptions which, when made, permit the inte- 
gration of equation (23) in terms of simple functions. The first is to assume that the 
cubic has two equal roots; the denominator then reduces to (z - b)\[z, - a and can be 
integrated. The two equal roots can occur only for specific values of o Q and q 0 , and 
for this reason, is of very limited usefulness. However, one of the cases of two equal 
roots gives the solutions for the A = 0 models. A second approach is to assume that 
density is zero. The assumption of zero density reduces the cubic to a quadratic which 
can be integrated directly. The solutions for the zero-density models are exact from a 
mathematical point of view, but, when used, it is assumed that the density of the universe 
is so small that the physical situation existing in the universe can be approximated by the 
zero-density model. The alternative approach to those already suggested seeks to obtain 
a mathematical approximation for equation (23). Such an approximation occurs when the 

3 

term 2ct q z in the cubic can be neglected with respect to the other terms in the equation. 

Methods for integrating equati on (23).- If equation (23) is written in the form 


o> = 


HqRq Jq 


dz 

\j(z - Of) (z - /3)(z -rj) 


it can be integrated by the methods given in standard texts on elliptic integrals. Such 
solutions, in general, provide only numerical answers once the roots are known, and thus 
imply that a Q and q Q are known. The alternative is to integrate equation (23) by 
numerical methods on a digital computer. In the problem of analysis of observational 
data, Oq and q Q are the unknown quantities, and expressions for co of the form 
co * f^a Q ,q 0 ,zj are more useful. 

Two equal roots .- There are two separate cases when there are two equal roots. 

In the first, A ^ 0 and in the second, A = 0. In the first case A * 0, the variable u> 
is given by 


<JL) = 


c 

H 0 R(W2ao 



dz 



* \L , 9c, ° ' 

3o- 0 J]j 6<Tq 


(24) 
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which for k = 1 integrates to 


[ 


M _ log \) 6g o z + 9ff o - (% + i) - 3cr o - fio + *) (% + !) + gj/300 - (g Q + 1) ^ 

\[6a 0 z + 9cr 0 - (q Q + 1) + l/3^3a Q - (q Q + 1) ^9a 0 - (q Q + l) - l/3^3a Q - (q Q + l) 


and for k = -1 


co = 2 


tan 


■ 1 j 6 °'o z + 9g o - (% + x ) 




+ 1 - 3or n 


- tan“ 


CO 

o q 

( q o + X ) 

P\j 

f % + 

l-3a 0 


(26) 


The necessary condition for these solutions to exist is given by 

3a o - (% + 1 ) 3 - 27(T o 2 ( cr o - %) = 0 


(27) 


These pairs of a Q and q Q that satisfy equation (27) are shown in figure 1 for 
-4.0 = q 0 = 6.0 and 0 ^ cr 0 S 2.0. The solutions for co correspond to models with 
k = ±1. The solutions for k = 0 models are points and occur when cr 0 = q 0 = 0.5 and 
a Q = 0, q Q = -1.0. Equations (25) and (26) do not apply at these points and special 
expressions are required for co. One of these points cr 0 = q Q = 0.5 belongs with the 
A = 0 models and is discussed later. The second point a Q = 0, q 0 = -1.0 is on the 
curve for negative q Q in figure 1. The co for this point is given by 


co = ■ 


HqR-O 


(28) 



Figure L- Combinations of the density and acceleration parameters, o 0 and q 0> that satisfy the 
conditions for the integration of equation (27). 
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In the second case A = 0, where equal roots exist, w is given by 


u> 


9 f 

H o R o f^ J 0 


dz 


(z + 1)1 z + 


%+ 1 


2 O n 


( 29 ) 


and the necessary condition for the solutions to be valid is that o Q = q Q ; by equation (8) 
this relationship corresponds to A = 0. In the case of A = 0 models, co is obtained 
in a more usable form if A is set to zero in equation (22) and then integrated as was 
done by Mattig (ref. 22). The resulting co equations are for k = 1, 


co = cos 


-1 


Ho^o^od + z ) 


- 1 


- cos 


-1 1 0“ 

, Hq^Ro^^o 


for k = 0, 


co = 


2c" 


2c" 


H 0 2 R 0 2 °b VVRqMI + z) 


2r% 2 , 


and, for k = -1, 


co = cosh 


-1/ c z 
Uo 2 Ro 2 »o 


+ 1 - cosh" 


, H 0 2 R 0 2 a 0 (l + z) 


+ 1 


(30) 


(31) 


(32) 


When A = 0, k = 0 only when q Q = 0.5. In A = 0 models, q Q must be positive or 
equal to zero since negative q Q would give negative density (that is, o Q < 0) ; and since 
objects (galaxies, quasi-stellar radio sources) exist in the universe, the density and con- 
sequently o 0 cannot be negative. In both cases of two equal roots A * 0 and A = 0, 
equation (9) is used to determine k and, consequently, determine the form of co 
to be used. 

Equations (25), (26), (28), and (30) to (32) represent exact solutions for co and 
impose a severe restraint on the relation between o Q and q Q . In both cases, o Q -A q Q 
and Oq = q Q , it is necessary for o Q and q Q to have a relationship that comes solely 
from the mathematics of the problem and is not a result of observation or the physics of 
the problem. 

The zero-de nsity model.- Observational results (ref. 16) indicate that the density 
of the universe is probably close to lO -2 ^ g/cm 2 . This level of density means that a 0 
is about 0.05, and for z < 1.0 the terms involving cr 0 will contribute little to the cubic 
in equation (23). It is therefore interesting to investigate the zero-density models, which 
are exact mathematical models, as approximations for a very low-density universe. 

When o Q is set equal to zero, equation (23) becomes 
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A 


Figure 2 compares the quantity coH 0 R 0 /c as calcu- 
lated by equations (23) and (33). The use of coH 0 R 0 /c 
permits the effect of the zero-density assumption on the 
integral to be determined. The differences between the 
exact and zero-density models are small to z = 1.0 

Hq 

for q Q = 0.5 and decrease with increasing q 0 . For ' 

z > 1.0, the differences increase. This increase occurs 
because the terms involving density are z^ and 7? 
and the contribution of these terms to the cubic in equa- 
tion (23) builds up very quickly above z = 1.0. The 
good agreement obtained between coH 0 R 0 /c for the 
exact and zero-density solutions indicates that for low Figure z- Differences in for the exact 

density, the zero-density model could be used in place fin ! ,e densit y and zero-density models of the 

of the exact model. When a simple model is used to 
replace a more complicated exact model, it is always 

good to have some idea of the differences between the simple and exact model. System- 
atic calculations with the exact and zero-density solutions showed that if the following 
inequality was satisfied 

<x 0 § 0.075q o (q o - l)+0.1 (34) 

the difference in coH 0 R 0 /c for the exact and zero-density solutions would be 2 percent 
or less at z = 1.0. 



Integration of equation (33) gives a general solution for co which is 


1 h q r 0 

CO = log 

fk 


(1 + z) + ll 

1-^(1 + s) 2 - q„ 

II 

B 

o 

0 

o 

t> 

c 

fE £.+ \ 

'S c2 - q 

HqRq 1 

„2 r 2 q o 
n 0 tt 0 


and S(co) is given by 


s( "> = + »)(» - % - o + 


(35) 


(36) 


As shown in figures 3 and 4, the agreement between the exact and zero-density 
solutions for co and S(co) is poor compared with the agreement for the corresponding 
solutions for coH Q R 0 /c. 
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Figure 3.- Differences in co for the exact finite density and zero-density 
models of the universe. (Eqs. (23) and (35).) 


Both co and S(co) (see eqs. (23), (35), and (36)) are multiplied by a coefficient 
that contains l/R 0 . In reference 1 it was shown that the assumption of zero density 
introduced large differences in R 0 . In addition to the differences in the integral, which 
were found to be small (see fig. 2), the only term that can introduce additional differ- 
ences, such as those shown in figures 3 and 4 as compared with figure 2, is R 0 . This 
relation suggests that if the zero-density result is multiplied by R 0 for the zero- 
density model and divided by R 0 for the exact model, the differences between the two 
models should approach the differences shown in figure 2. Accordingly, co* and S*(co) 
to replace co and S(co), respectively, of the zero-density model are defined as follows: 


CO* = CO. 


Wd 

1 ( R °)e 


(37) 




( 38 ) 


16 


2.2 r 



Figure 4.- Differences in S(w) for the exact finite density and zero-density 
models of the universe. (Eqs. (21a), (23), and (36).) 

where the subscripts d and e stand for the zero density and exact models, respec- 
tively. As shown in figures 5 and 6, the differences between co* and a> and between 
S and S(cu), where o> and S(a>) are for the exact model, approach the differ- 
ences found for u>H 0 Ro/ c * 

The use of and S*(co) is required only when it is desired to investigate the 

variables by themselves or when c o and S(o>) are substituted in an equation that does 
not cancel the l/R 0 term that is introducing the large differences. Therefore, in equa- 
tions (14) and (15), S*(co) should not be used, whereas in equation (16) and S*(cu) 
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Figure 5.- Differences in w or uf for the exact finite 
density and modified zero-density models of the universe. 
(Eqs. (23) and (37).) 



z 


Figure 6.- Differences in SM or S*(to) tor the exact finite 
density and modified zero-density models of the universe. 
(Eqs. (21a), (23), and (38).) 


should be used. Lastly, it should be noted that the use of co* and S*(w) requires that 
an estimate of the average density of the universe be available. 

If the zero-density model is to be used for the analysis of observational data, the 
k and A associated with the zero-density model must have the same signs as the k 
and A of the nonzero-density model as long as o Q = q Q . This requirement means that 
zero-density models cannot be used for the analysis of observational data unless the fol- 
lowing conditions in the observed universe are present: For q Q > 0, the zero-density 

model can be used if q Q > a Q and q Q + 1 > 3a Q ; for -1 < q Q < 0, if 1 > 3cr 0 - q Q ; for 
q Q < -1, if there is no restriction on q Q . 

Approximate model so lut ions for co . The results shown in figure 3 for coH 0 R 0 /c 
indicate that the integral in equation (23) is very little affected by cr 0 . Therefore, an 
investigation was made of the effect of assuming that for small cr 0 and zSl, the term 

O 

2c q z° makes a negligible contribution to the cubic in equation (23). This assumption is 
attractive because the only differences introduced are in the integral wH 0 R 0 /c, because 
R 0 is not affected, and because the resulting solutions permit the determination of a 0 
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o 

from observational data. When it is assumed that 2o Q z is a term that makes a neg- 
ligible contribution to the cubic in equation (23), equation (23) becomes 



Figure 7 compares the wH 0 R 0 /c obtained from equation (23) and from equation (39). 
As can be seen, the difference between the solutions is practically nonexistent out to 
z = 1. 



0 4 .8 1.2 1.6 2.0 2.4 2.8 


H 

Figure 7.- Differences in — g-® for the exact and approximate finite density models 
of the universe. (Eqs. (23) and (39).) 


The integration of equation (39) for co gives 
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and substitution into equation (21a) for S(u>) gives 



(41) 


Figures 8 and 9 compare the co and S(u>) for the exact and approximate models. 
The exact solutions for co and S(u>) are given by equations (23) and (21a), and the 
approximate solutions by equations (40) and (41). The differences in co and S(co) for 
the exact and approximate models are small, on the order of those found for 0 >H o R o / c. 

As indicated in figures 7 to 9, this approximate solution for o> and S(co) has 
been checked for a Q = 0.08, and q 0 = 0.5 and 2.5. Charts are given in appendix A that 
show the o 0 and q Q combinations that give a 2- and 4-percent error in this approxi- 
mation for co. 



Figure 8.- Differences in w for the exact and approximate finite density 
models of the universe. (Eqs. (23) and (40).) 
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Figure 9.- Differences in S(w) for the exact and approximate finite 
density models of the universe. (Eqs. (21a), (23), and (41).) 


Another method to obtain an approximation for oo is to expand the integral in 
a McLaurin series about z = 0 and then expand the expression for S(o>) in an appro- 
priate series. This form of the red-shift — magnitude relation, called the expansion 
form, has been extensively studied (see ref. 15), and the red-shift — magnitude relation 
using this approach is 



This relation is restricted by assumptions made in its development to small values of z, 
that is, z « 1. Sherman (ref. 1) and Mattig (ref. 22) have compared this form of the 
red-shift — magnitude relation with the zero-density and A = 0 red-shift — magnitude 
relations and found differences of 1.8 magnitude and 0.4 magnitude, respectively. As 
the zero-density model and the A = 0 model are special solutions of the exact case, 
there is no reason to expect better agreement between the exact and the expansion model 
than has been demonstrated in the special cases. 


The Scale Factor and Time Equations 

There is other information that can be obtained from the model universe. This 
information is the scale factor, time of light travel, and time since the beginning of 
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expansion. Equation (7b) is integrated to obtain either the time of light travel or the 
time since the beginning of expansion. 


The starting point for these calculations is equation (10) which is equation (7b) in 
integral form. Equation (10) in the form 


t 0 - t 


'“o' 1 ! 


dz 


0 (1 + z) 2a 0 z 3 + (3cr 0 + q Q + l)z 2 + 2(q 0 + l)z + l] 


i/2 


(43) 


is used to determine the time of light travel and in the form 


T = 





dy 

+ (l + % 



(44) 


is used to determine the time since the beginning of expansion. In equation (44), T is 

the time since the beginning of expansion and y = — . 

K 0 


Time of Light Travel 


Equation (43) is integrated to obtain the time of light travel. The solution to equa- 
tion (43) is of the form t Q - t = f(R). This solution is then inverted to give the scale 
factor R = h (t Q - t) . 

The case of two equal roots .- The cubic in equation (43) has equal roots for certain 
combinations of a Q and q Q when A * 0 and always when A = 0. When A * 0 and 
the cubic has two equal roots, equation (43) becomes 


^0 



dz 


(1 + z)( z + 


q ° + ~ ( q p ~ s 

3cr 0 /| 6cr 0 


(45) 


This integral may be written in the form 

-1 


t = 


3 °bH 0 


-r 

+ 1 - 3CT o) J ° 



1 



1 

(z + 1 )| 1 

1 9 a o 

1 ry I 


L + q ° + \ 

L , 9<I o - (q„ + 1) 

6°b 

{ 3a 0 )\ 

| 6a 0 


dz 


(46) 


and in this form the equation for t 0 - t can be easily integrated by standard formulas. 
The integrals in equation (45) or (46) give t 0 - t only when o 0 and q 0 satisfy the 

— O O 

condition j3a Q - ^q Q + 1)J - 27cr 0 (a 0 - q Q ) = 0. Figure 2 shows a plot of the values of 
q 0 and a Q that satisfy this condition. In the case of t Q - t, the expression resulting 
from the integration of equation (46) is not as simple as it was for co and, in general, 
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will consist of two inverse tangent functions, or two logarithmic functions, or a combina- 
tion of inverse tangent and logarithmic functions. There are 11 possible combinations of 
these functions and the specific combination depends on the o Q and q Q used in the 
integration of equation (46). Because of the many different conditions that would have to 
be covered and the highly restricted nature of the solutions, the integration of equa- 
tion (46) is not presented in detail. However, it is easily integrated by using the stan- 
dard formulas. The second case of two equal roots occurs when A = 0. When A = 0, 
two roots are equal and equation (43) becomes 


t = 


H „ _1 rz 




Jr 


dz 


0 (1 + z) 2 i/z + 


1 

2q 0 


(47) 


and the necessary condition is that a Q = q Q . In this case it is simple to use equation (10) 
with A = 0 to obtain 

to-t=f Ro - ^ dR (48) 

J R 


(2H 0 2 R 0 3 q 0 - kc 2 Rj 


2 ^ 1/2 


Because of difficulties that arise when inverting the solution of equation (48) to obtain 
the scale factor, it is convenient to introduce the dummy variable ? when integrating 
for k = -1 or 1. When k = 0, the inversion of the integration of equation (48) presents 
no difficulties. When k = -1, the substitution c 


R = sinh 3 ? is used and 


2H 0 2 R 0 3 q 0 

R = sin 2 ? is used. The results of the inte- 


when k = 1, the substitution — 

2H 0 2 R 0 3 q 0 

gration of equation (48) for k = -1 and k = 1 by using the substitutions, and the 
direct integration when k = 0 are as follows: 


^ " t = 


^ - 7 -r c osh(? 0 + ?)sinh(?o -?)+?- ?o) (k = -1) 

H 0 (l - 2q 0 f /2L J 


R = 


«l 


o( i.i o) »^ coah2C ‘ 1) 


> (49) 


(k= -1) 


to_t= fer 


i - jif 72 


R = R 0 \l - 


r 

2/3 



2 (‘o’ 1 ) 

J 


(k = 0) 


(k = 0) 


/ (50) 
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to - t = — - s + sin (?o - ?)cos(? 0 + d (k = 1) 

H 0 (l - 2q 0 ) 3 / 2 L J 


> (51) 


R = 


c % 


H 0 (2q 0 - 1) 


372 


(1 - cos 2£) 


(k = 1) 

j 


The foregoing equations for t 0 - t and R are not too useful in that a knowledge 
of one of the variables is required in order to determine the other. In the case of R 
for k = ±1, the equations are transcendental as £ is a function of R. However, all the 
foregoing equations may be expressed as a function of the red shift and acceleration 
parameters through the use of the relation = 1 + z and equation (8) . When using 

rv 

equation (8) for the A = 0 case, o Q is set equal to q 0 . 

The zero-density m odel.- In the case of the zero-density model universe, equa- 
tion (10) becomes 


to - t = 





(52) 


The integration of equation (52) gives the general expression for t Q - t for the zero- 
density model which is 

H„(t„ - t) = Jq^cos- 1 |L0I - eo,^) (53) 

Figure 10 compares t Q - t for the zero-density model universe with the t 0 - t for the 
general case which is given by equation (43). As can be seen, there is excellent agree- 
ment between the two solutions where cr 0 = 0.08 and 0 and q Q = 0.5 and 2.5. 

The scale factor R is obtained by inverting equation (53) so that R is given as 
a function of t 0 - t. This inversion gives for the zero-density model 


R — Rr 


cos 


ll^o-t) 


H, 


-1 




sin 


\M fc o ■ *) 


H r 


(54a) 


and is a general expression for R for the zero-density model universe. Figure 11 
compares solutions for R for the exact and zero- density models and the agreement is 
poor. If a modified zero-density scale factor R* is adopted 


R* = 


a- 


(54b) 
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Figure 10.- Differences in H 0 (t 0 - t) for the exact finite density and zero-density 
models of the universe. (Eqs. (43) and (53).) 


as was done for w and S(o>), as shown in figure 12, there is very good agreement 
between the exact and zero-density modified scale factors. 

The zero-density model universe and the A = 0 model universe have a common 
model and that model is for q Q = 0. The restrictions on the use of oo* in the zero- 
density model apply to the use of R* when it is used to determine t Q - t and R. 

Approximate solutions for t Q - t and R.- The approximate solutions for t Q - t 

and R make use of the same approximation that was used for obtaining the approxi- 
mate solutions for co, that is, the term 2a Q z^ makes a negligible contribution to the 
integral in equation (43) . When the term 2a 0 z^ is assumed to make a negligible con- 
tribution to the cubic in equation (43), equation (43) becomes 


t 


o 



dz 

(1 + z)f(3cr 0 + q Q + l)z 2 + 2(q 0 + ljz + 



(55) 
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Figure 11.- Differences in H 0 R for the exact finite density and zero-density 
models of the universe. (Eqs. (43) and (54).) 





and integrates to 


Ho " 1 

. -l ( 3 cr o-%) R - 3 a o R o .-1 

sin - sin 

\j% ~ 3o o 

R oK( q o + x ) - 3a o fio(% + - 3o o 


( 56 ) 


As shown by figure 13, equation (56) is an excellent approximation for t Q - t given by 
the exact expression (eq. (43)). Inversion of equation (56) gives the expression for R 
which is 


R = 


R 0 


sini 

fq ° " 3a °/ t - t ) 

3a 0 R 0 

\/% ~ 3a o 

_K- 3a o h 0 -i 


Ho' 1 ( ° 

% - 3a O 


(57) 


As in the case of o>, t Q - t and R obtained from this approximation are not 
applicable to the A = 0 models. For negative values of q Q , the results obtained from 
the approximate solution are not as good as those for q Q > 0. The approximation 
becomes worse when q Q = -(l + 3cr 0 ), and the approximation improves when q Q 
becomes more negative. 



Figure 13.- Differences in H 0 (t 0 - t) for the exact and approximate finite density 
models of the universe. (Eqs. (56) and (43).) 
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Charts that show these combinations of a Q and q 0 that give a difference between 
the exact and approximate solutions of 2 percent and less or 4 percent and less are pre- 
sented in appendix A. 


Time Since the Beginning of Expansion 

Equation (44) gives the time since the beginning of expansion for those models of 
the universe that have R = 0 at t = 0 and excludes models 10, 11, and 14 of table I. 

If the proper initial value of R at t = 0 is used for the lower limit, a time since the 
beginning of expansion can be calculated for models 11 and 14. In the case of model 10, 
time since the beginning of expansion is meaningless because this model is a static 
model and there is no expansion. 

Equation (44) is a definite integral, and it can be easily integrated for the zero- 
density model and the A = 0 model, the significant family of models when two roots are 
equal. In the case of the approximate model, equation (44) takes the form 




dy 


( 3ct o - * o ) y 2 - 6a 0 y + ( 3ct o + % + x ) 


n 1 / 2 


(58) 


which can also be easily integrated in terms of simple functions. 

Calculations of the time since the beginning of expansion are mainly useful in deter- 
mining whether the time scale inequality which is 

Time since beginning of expansion i Age of oldest stars (59) 


is satisfied for a given model of the universe. Equation (59) is meaningful only when a 
model starts from the point R = t = 0. At this point it is assumed that no stars and 
galaxies existed and that they were formed shortly after expansion started. Thus, time 
since the beginning of expansion has significance in these models. However, for 
models 11 and 14 that have open negative time axes and finite scale factors, equation (44) 
loses its meaning because stars and galaxies could have existed before expansion started. 


DISCUSSION 


Four simple models of the universe have been derived and equations that relate 
observation and theory have been obtained for each model. Of the four simple models, 
only two are significant because fewer restrictions are incorporated. These models are 
the zero-density and approximate model. Both of these models allow all values of k 
and A and if the density is on the order of that predicted by observation or less 
(p ~ 10 _3 0 g/cm 3 or less), then these models are highly useful. The other two models, 
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those that have two equal roots, are not considered significant because of the restric- 
tions placed on a Q and q Q by the condition of two equal roots. However, one of these 
two models, the model with A = 0, has been used in the past (ref. 19) and therefore is 
compared with the zero-density model. One other model, the expansion model (eq. (42)), 
is considered in the comparison because of its frequent use. 
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Figure 14.- Red-shift— apparent-magnitude data used for the 
comparison of models. 


It was decided to use the red- 
shift — magnitude relation for the exact 
zero density, A = 0, approximate and 
expansion models to determine <J 0 and 
q Q from a set of observational data. 
Two sets of observational data are 
available: one by Humason, Mayall, 

and Sandage for 18 clusters of galaxies 
and one by Baum for 8 clusters of gal- 
axies. Both of these sets of data are 
summarized in reference 3. Because 
Baum's data (fig. 14) contain six color 
photoelectric magnitudes, no K cor- 
rection is required. For this reason 
and because Baum's data appeared to 
have less scatter, it was decided to use 
it for the calculations to compare the 
various models of the universe pre- 
sented in this paper. The results of the 
calculations are given in table II. 


TABLE H.- RESULTS OF CALCULATIONS 


Model used 

Equation 

CT o 

% 

k 

A 

Model (Table I) 

Exact 

(14), (21a), and (23) 

3.029 

0.0077 

1 

>0 

9 (Le Maitre) 

Approximate 

(14) and (41) 

3.20 

-.0059 

1 

>0 

9 (Le Maitre) 

A = 0 

(14) a =Q , (21a), (30), (31), and (32) 
o ^0 

1.56 

1.56 

1 

0 

8 (Oscillating) 

Zero density 

( 14 ) cr =0 and (36) 

0 

2.934 

-1 

<0 

1 (Oscillating) 

Expansion 

(42) 

— 

1.39 

-1* 

<0* 

1 (Oscillating) 


^Equations (8) and (9) were used to compute k, and A for this model at a value 
of 0.45 (p Q ~ 10 - 30 g/cm^j was used. 
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The results presented in table II are the initial results of a study of the problem of 
determining cr Q and q Q from observational data and are results for a very limited 
amount of data. The model characteristics that are indicated when density is neglected 
(expansion model) or a restrictive assumption made (A = 0 and zero-density models) 
disagree with the model characteristics that are indicated when the exact red- shift — 
magnitude relation is used. However, the ct 0 and q Q predicted by the approximate 
model are close to those predicted by the exact model. This result is not surprising 
since the cr 0 and q Q obtained through the use of the exact model lie in a region where 
the approximate model is a good approximation for the exact model if an error between 
2 and 4 percent in the elliptic integral is tolerable. Both the exact and approximate 
models predict a 0 on the order of 3.0; this prediction represents a density that is about 
20 times higher than the maximum density based on observation obtained by Oort (ref. 13) 
and about 60 times Abell's best estimate (ref. 12). If the trends established by these 
initial results are substantiated by further calculations with more extensive data, only 
the exact and approximate models should be used for the analysis of observational data. 
Solheim's results (ref. 23) appear to confirm these results. 

The results presented in table II have been plotted in figure 15 with the data super- 
imposed on the curves. The most striking fact about this figure is that for the range of 
z covered by data, 0.004 to 0.460 (that is, log z = -2.398 to log z = 0.361), all the 



m c 

Figure 15.- Comparison of several model universes with red shift z and apparent 
magnitude m c for Baum's eight clusters from reference 3. 
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models of table II appear to fit the data equally well. Above z = 0.5 (log z = -0.30103), 
the various models diverge. Present data for galaxies do not extend above z = 0.5, and 
data with z >0.5, based on the evidence of figure 15, are required to determine a model 
of the universe. 

Quasi- stellar radio sources (see ref. 24) have measured red shifts greater than 
0.5; however, these sources were not used to extend the data for galaxies used in the 
calculations because (1) the cause of the red shift is in doubt and (2) the intrinsic lumi- 
nosity of these sources appears to vary widely and is different from that of the galaxies. 
The use of the red- shift — magnitude relation requires the assumption that all sources 
have about the same intrinsic luminosity. Additional data from quasi- stellar radio 
sources published since reference 24 have not clarified the situation. 

CONCLUDING REMARKS 

Under the assumption that general relativity is the theory of gravitation that best 
represents the observed universe, relativistic models that are both isotropic and homo- 
geneous have been studied. In addition to the assumptions of isotropy and homogeneity, 
it was necessary to assume that pressure terms could be neglected with respect to den- 
sity terms in order to solve the field equations that describe the universe. 

The solutions for the scale factor and radial metric variable to, which are gen- 
erally elliptic integrals, have been studied for the exact case and several simplified 
cases. In one of the simplified cases, the integrals for the scale factor and radial metric 
variable have two equal roots and although both integrals can be integrated by simple 
functions, the condition of two equal roots imposes too severe a restriction on the pos- 
sible values of the density and acceleration parameters for the solutions with two equal 
roots to be of more than academic interest. Models with the cosmical constant equal 
to zero are contained in the case of two equal roots and the same remarks apply. In 
addition, the cosmical constant is a constant of integration in the field equations and 
should not be arbitrarily set to zero. Simplified models that have two equal roots should 
not be used unless the analysis of observational data with a more general model indicates 
that the special conditions for the integrals for the scale factor and radial metric variable 
to have two equal roots exist. 

Zero-density model universes and an approximate model in which the cubic term 
was assumed to make a negligible contribution to the integral in the definition of the 
radial metric variable co were also studied to determine their applicability to the anal- 
ysis of observational data. 

A computing program based on the method of differential corrections and the data 
for Baum's eight clusters of galaxies was used to evaluate the exact, approximate, A = 0, 
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zero-density, and expansion models. Both the exact and approximate models predict a 
high-density universe and the acceleration parameter for these models is not in agree- 
ment with that determined by the A = 0, zero-density, and expansion models. Because 
of the high density and the lack of agreement between the values of the acceleration 
parameter determined by the various models, except in the case of the exact and approxi- 
mate models which predict an almost zero value for the acceleration parameter, it was 
concluded that only the exact and approximate models should be used for the analysis 
of observational data. The approximate model should be used only if an error between 
2 and 4 percent in the integral in the radial metric variable co is considered tolerable. 

Further study of the results obtained from the computing program showed that in 
the range of red shift covered by present data, all the models used appeared to fit the 
data equally well, and data that include red shifts greater than one-half are required to 
define a model of the observed universe. 

Although the data for quasi- stellar radio sources fall in the required range of red 
shift, it was pointed out that, as yet, these data are not considered in the problem of 
determining a model of the universe because of the uncertainty as to the cause of the red 
shift and because the intrinsic luminosity of quasi- stellar radio sources shows a wide 
variation and is different from that of galaxies. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., September 21, 1967, 

188-41-01-01-23. 
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APPENDIX A 


DETERMINATION OF REGIONS OF APPLICABILITY 
OF THE APPROXIMATE MODEL 


The region of CT 0 ,q 0 over which the approximate model is applicable is deter- 
mined for errors of 2 and 4 percent introduced through use of the approximate model. 

The only difference between the approximate and exact models is the approximation of 
equation (23) by equation (39) for w and the approximation of equation (48) by equa- 
tion (55). These equations were set up as inequalities, and a high-speed digital computer 
was used to determine ct 0 and q 0 as a function of z. The following inequalities were 
used: 


K z ) 3 9 - I( Z ) 2 3 ^ 

!(z) 2 3 


(Al) 


for (jo and 

I ^55 " ^ z) 48 < 
«z>48 = 6 


(A2) 


for t Q - t. In these equations, I(z) indicates the integral and the subscript denotes the 
equations from which it came. Two values of the error e, 0.02 and 0.04, were used. 


The results of the calculations are presented in figure 16 for to and in figure 17 
for t Q - t. When e was 4 percent, it was found that the error for z = 0.5 was always 

less than 4 percent for the range of a Q and q Q considered; thus, no boundary for 

z = 0.5 is shown in figures 16 and 17, the region between the axis and the lines of con- 
stant z is usable, and a combination of o Q and q Q that falls on or inside of a line of 

constant z is usable. For instance, if ct q and q Q fall on the z = 1.0 line, the 

approximate model could be used to determine cr 0 and q Q out to z = 1.0. The maxi- 
mum error would occur on the z = 1.0 line; if a Q and q 0 fall inside this line, the 
error is less. 


In practice, the charts would be used as follows. The approximate model would be 
used to determine cr Q and q Q , and then the determined o Q and q Q and the maximum 
z used would be entered into the chart to determine whether the approximation is valid. 
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APPENDIX A 



(a) £ = 0.02. 


.0 1.5 



i i _!_ .1 I 1 
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Ob 


(b) £ = 0.04. 

Figure 17.- Regions in which the approximate model 
gives valid solutions for tg - t. The useful region 
is between the q 0 axis and lines of constant z. 
(Eq. (A2).) 



APPENDIX B 
SYMBOLS 

a = 7.564 x KT 15 

a,b roots for cubic with two equal roots 

C = 5 log 1Q + M Q - 5 
n o 

c speed of light in vacuo 

D c distance by measuring rod, R 0 S(w) 

Dj luminosity distance 

E curvature scalar 

E^y contracted Riemann-Christoffel tensor 

f( ) function 

G universal constant of gravitation 

g^j, metrical tensor 

H 0 Hubble parameter 

h inverted function 

I(z) integral 

K red- shift correction 

k curvature constant 

M 0 absolute magnitude of equivalent local source 

m apparent magnitude 

m c corrected apparent magnitude 
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APPENDIX B 


N(m) 

n 

P 

Q 


% 

R 

Rl,R2 

Re 

R* 


S(w) 

S*(co) 


T 

t 

tn 


v 

V ds 


V» = 


dx^ 

ds 


number of like sources brighter than a given magnitude 
number of like sources per unit volume 
total pressure 

number of square degrees in celestial sphere 
acceleration parameter 

scale factor that describes geometric history of universe 

roots of cubic in equation (10) 

value of R associated with Ae 

value of R for modified zero-density model 

coordinate in metric subspace 

function of co that depends on curvature of space 

value of S(u>) for modified zero-density model 

time along line element 

radiation temperature, °K; also time since beginning of expansion 
time 

present time 


v 0 random radial velocity 

xP,x v coordinates 


y 


R_ 

R 0 
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APPENDIX B 


z 

a,P,V 


C 


6 



A 

A e 

I 

P 

P o 
°o 

a) 

co* 


red shift 

dummy roots of cubic 
error 

linear diameter; also dummy variable 
coordinate in metric subspace 

cosmical constant 

critical value of cosmical constant 

function of radial metric variable r 

material density 

present material density 

density parameter 

angular diameter; also coordinate in metric subspace 

function of radial metric variable 

value of co for modified zero-density model 


Subscripts: 


e,d exact and zero-density models 

o present values 

/!, v tensor indices (0, 1, 2, and 3) 


Dots over symbols denote differentiation with respect to time. 
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